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Abstract. 

The study of photoionised gas in planetary nebulae (PNe) has played a major role in the achievement, over the years, of a 
better understanding of a number of physical processes, pertinent to a broader range of fields than that of PNe studies, spanning 
from atomic physics to stellar evolution theories. Whilst empirical techniques are routinely employed for the analysis of the 
emission line spectra of these objects, the accurate interpretation of the observational data often requires the solution of a set of 
coupled equations, via the application of a photoionisation/plasma code. A number of large-scale codes have been developed 
since the late sixties, using various analytical or statistical techniques for the transfer of continuum radiation, mainly under the 
assumption of spherical symmetry and a few in 3D. These codes have been proved to be powerful and in many cases essential 
tools, but a clear idea of the underlying physical processes and assumptions is necessary in order to avoid reaching misleading 
conclusions. 

The development of the codes over the years has been driven by the observational constraints available, but also 
compromised by the available computer power. Modern codes are faster and more flexible, with the ultimate goal being 
the achievement of a description of the observations relying on the smallest number of parameters possible. In this light recent 
developments have been focused on the inclusion of new available atomic data, the inclusion of a realistic treatment for dust 
grains mixed in the ionised and photon dominated regions (PDRs) and the expansion of some codes to PDRs with the inclusion 
of chemical reaction networks. Furthermore the last few years have seen the development of fully 3D photoionisation codes 
based on the Monte Carlo method. 

A brief review of the field of photoionisation today is given here, with emphasis on the recent developments, including the 
expansion of the models to the 3D domain. Attention is given to the identification of new available observational constraints 
and how these can used to extract useful information from realistic models. 

Keywords: <Enter Keywords here> 

PACS: <Replace this text with PACS numbers; choose from this list: http : / /www . aip . org/pacs/ index . html> 



INTRODUCTION 

Photoionised plasma is present in many astrophysical en- 
vironments, from H II regions and Planetary Nebulae 
(PNe), that mark the beginning and the end stages of stel- 
lar evolution, to ionised interstellar and intergalactic me- 
dia, to the gas photoionised by high energy sources in 
AGNs and seen in distant quasars. The interpretation of 
their complex emission line spectra relies on our ability 
to disentangle the underlying physics, which involves a 
number of microscopic atomic processes, highly sensi- 
tive to the physical properties of the emitting gas, and to 
the radiation field of the ionising source(s). Whilst em- 
pirical studies of the observations are routinely carried 
out to extract some basic information from the spectra, 
the application of large-scale numerical codes is often 
essential to the understanding of these sources. 

PNe, as many photoionised gas clouds, are generally 
transparent to the emission line radiation that is pro- 
duced in their interior and, therefore, a prediction of 
their spectra simply requires a knowledge of the elec- 
tron temperature and ionisation structure at each position 
in the ionised volume. This can be achieved by numer- 
ically solving the coupled ionisation balance equations 
and the thermal equilibrium. The local radiation field, 



including the stellar and diffuse components, must first 
be estimated by, classically, analytical or, more recently, 
statistical means. All major heating and cooling chan- 
nels must be accounted for as well as all ionisation and 
recombination processes. A comprehensive reference is 
the textbook by Osterbrock (1989), additionally, a review 
on quantitative spectroscopy of photoionised clouds was 
given by Ferland (2003), whilst a summary of emission 
line analysis can be found in a number of papers pre- 
sented at a conference to honour Silvia Torres-Peimbert 
and Manuel Peimbert on their sixtieth birthdays (Henney 
et al. 2002) and those presented at a previous meeting 
organised to celebrate the seventieth birthday of Don Os- 
terbrock and Mike Seaton (Williams & Livio 1995). 

The first generation of codes that were able to treat 
the temperature stratification in H II regions appeared 
in the late 1960s (e.g. Flower 1968, Rubin 1968, Har- 
rington 1968), these codes were severely hindered by the 
limited computing power available at the time and the 
lack of a comprehensive atomic data set. The fast im- 
provements in hardware and the large volume of atomic 
data become available in the last forty years have allowed 
photoionisation calculations to become more and more 
sophisticated boosting the predicting power of the codes 
and making them an essential tool for modern quanti- 



tative spectroscopic studies. At present more than eight 
large scale photoionisation codes are in use. Their perfor- 
mance has been tested through a number of workshops, 
resulting in the definition of a set of rigorous bench- 
marks (Pequignot et al. 2001), generally referred to as 
the Meudon/Lexington benchmarks. The state of the art 
of the field is accurately described by a number of papers 
presented at the latest of these meetings and included in 
Spectroscopic Challenges of Photoionised Plasmas (Fer- 
land & Savin 2001). In the present paper I summarise 
where photoionisation modelling is today, highlighting 
the major developments that have occurred in the past 
five years, with particular emphasis on the applications 
of modern plasma codes to the modelling of PNe. 

NEW ADVANCES: HOW? 

The construction of a photoionisation model aims to 
match all available constraints for a given target, which 
may include spatially resolved or not, optical, infrared 
(IR) and/or ultraviolet (UV) spectra. Images in one or 
more narrow-band filter are also available for many ob- 
jects. A well-constrained tailored photoionisation model 
can unveil a large volume of information regarding the 
nebular properties, including the gas density distribution 
and its chemical composition, and regarding the proper- 
ties of the ionising star(s), such as effective temperature 
and surface gravity. Additionally, models of PNe have 
been used in the past to constrain the magnitude of some 
then unknown atomic process rates, as in the case of the 
hydrogen-oxygen charge exchange process (Pequignot, 
Stasinska & Aldrovandi, 1978). 

The desire of obtaining realistic models for a variety 
of emission line sources is the main drive behind the 
continuous developments of large photoionisation codes 
in the last forty years, including the more recent advances 
that will be discussed here. Progress has been made in 
many aspects, in this paper the follow five main streams 
will be considered: atomic data updates, the treatment 
of time-dependent effects, the inclusion of dust grains, 
the expansion of the models to the PDR and, finally, the 
development of 3-D codes. These will be discussed in 
more detail in the remainder of this section. 



Atomic data updates 

Advances in computer hardware and the construction 
of large scale theoretical atomic codes that can compute 
photoionisation cross-sections, electron-ion recombina- 
tion rates, level energies, electron-collision strengths, 
transition probabilities and oscillator strengths for a large 
number of ions and transitions have allowed a compre- 



hensive data set to be compiled. In particular, the Opac- 
ity Project and subsequently the Iron Project (Berrington 
et al. 1987), have produced a wealth of astrophysically 
useful atomic data, that is readily accessible via the TIP- 
TOPbase website (see e.g. Mendoza et al. 2002, Nahar 
2003). 

New total recombination coefficients, which include 
both radiative and dielectronic contributions, have re- 
cently been presented for a number of ions (e.g. Nahar 
2005, and references therein). These are ab-initio cal- 
culations that simultaneously and self-consistently treat 
photoionisation and recombination processes, although 
in some cases they rely on theoretical predictions of 
the resonance positions. The emphasis has recently been 
on higher ionisation species, not expected in PNe. New 
calculations of total recombination coefficients of some 
third-row elements in ionisation stages commonly ob- 
served in PNe are currently under way (Storey, Ercolano 
& Badnell, in preparation). These use the AUTOSTRUC- 
TURE program of Badnell (see e.g. Badnell 1999) in the 
LS- or intermediate-coupling (IC) schemes, where nec- 
essary, with corrections to the theoretical level energies 
being applied, where experimental values are known. 
The coefficients are calculated for a wide temperature 
range and can also be applied to very cool (approx. 
500 K) plasma. 

Observations of optical recombination lines (ORL) in 
long-slit and integral field unit (IFU) spectra are nowa- 
days readily obtainable. The analysis of such spectra 
poses some of the most intriguing questions yet to be an- 
swered in nebular astrophysics, namely the fact that ionic 
abundances derived from ORLs are consistently higher 
than those derived from the CELs in the same spectra, 
whilst the former diagnose temperatures that are signifi- 
cantly lower than those estimated using CEL ratios, He I 
recombination line ratios or the Balmer decrement (see 
e.g. Rubin 1997, Liu et al. 2000, Liu et al. 2001b, Tsamis 
et al. 2004). Photoionisation models of clouds contain- 
ing chemical inhomogeneities and/or density discontinu- 
ities may be crucial to clarify this long-standing problem 
(e.g. Pequignot et al. 2003, Tylenda 2003, Ercolano et al. 
2003), however they rely heavily on accurate recombi- 
nation data extending to low-temperatures. Work is cur- 
rently being carried out at UCL for a number of ions and 
some preliminary results have been presented by Bastin 
& Storey (2005a,b these proceedings). 



Time-dependence effects 

The ionisation structure and electron temperature of 
a photoionised cloud is obtained by simultaneously 
solving the thermal balance and ionisation equilibrium 
equations, classically under the steady-state assumption. 
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FIGURE 1. (van Hoof et al. 2004) The variation of the line ratio to H/3 for selected emission lines as a function of the grain size 
distribution. The left panel shows the results for the H II-region models, the middle panel for the PN models with silicate, and the 
right panel for the PN models with graphite. All line flux ratios are normalised to the base model without dust. The models numbers 
are defined in Table 4 of van Hoof et al. (2004). 



This approximation is valid when the atomic physics 
timescales are short compared to those of gas-motion or 
to the rate of change of the ionising field. If this is true 
the photoionisation problem can be considered to be fully 
decoupled from the dynamics; in this case a given gas 
density distribution represents a snapshot of the cloud 
evolution and a time-independent photoionisation sim- 
ulation can be performed. PNe and H II can generally 
be expected to be mostly in equilibrium, with dynamical 
effects only becoming significant in a very small region 
near the ionisation front, hence having near negligible ef- 
fects on the integrated emission line spectrum (e.g. Har- 
rington, 1977). 

The cases studied by Harrington (1977) were lim- 
ited to weak D-fronts with subsonic gas motions with 
respect to the front. There are, however, a number of 
instances when dynamics may have substantial effects 
on the predicted spectrum. Transonic photo-evaporation 
flows, such as those believed to be occurring in cometary 
knots of PNe (Lopez-Martin et al. 2001), cannot be ac- 
counted for under the assumption of thermal and ionisa- 
tion equilibrium. 

In the cases when shocks are present the gas cannot 
be considered to be in equilibrium, and classical time- 
steady codes are not appropriate. The MAPPINGS III 
code (Sutherland & Dopita 1993) was designed to deal 
with shock ionisation and has been applied to date to 
a variety of astrophysical environments. This code have 
recently undergone substantial updating and it is publicly 
available on-line via an interactive web interface (see 
http://cfa-www.harvard.eduMkewley/Mappings/i. 

Time-steady dynamics have recently been added to the 
widely used plasma code CLOUDY (publicly available 
at http://www.nublado.0rg/1; the numerical algorithms 
employed and a set of self-consistent dynamic models 
of steady ionisation fronts are presented in recent papers 



by Henney et al., 2005 a,b. 



Treatment of dust grains 

The importance of the effects of dust grains that are 
mixed with the gas in a photoionised region has long 
been known (e.g. Spitzer 1948), with more recent stud- 
ies also aimed at investigating the role of dust grains in 
the thermal balance of PNe (e.g. Borkowski & Harring- 
ton 1991, Ercolano et al. 2003c, van Hoof et al. 2004). 
The grains are heated via the absorption of UV photons 
from the continuum as well as resonance emission line 
photons, for example H I Lya, C iv A 1549, N v A 1240, 
C II A 1336, Si IV A 1397, Mg II A2800. The absorbed ra- 
diation is mainly re-emitted by the grains in the IR, and 
this provides the main dust cooling mechanism. Addi- 
tionally, the gas and dust components are coupled by a 
host of microphysical processes, including photoelectric 
emission from dust grains (for a recent discussion see 
Weingartner & Draine 2001), which may be an impor- 
tant gas-heating mechanism in H-deficient environments 
and in PDRs. Gas-grain collisions provide an extra cool- 
ing channel for the gas and a heating channel for the dust, 
and ionic recombination on dust grains may also have a 
small effect on the ionisation structure of the gas. 

Over the last thirty years, a number of photoionisation 
codes have included some kind of treatment to allow for 
dust opacities to be taken into account for in the ionised 
region and/or obtain estimates for the dust photoelectric 
heating. This topic, however, has recently been revisited 
by some authors in view of the expansion of existing 
codes to allow the construction of self-consistent models 
of ionised regions with their surrounding PDRs, where 
dust provides the main source of opacity. This will be be 
discussed in more detail in the next section. An exam- 
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FIGURE 2. Ionisation structure of a standard cloud from Abel et al. 2005; the two vertical lines represent the hydrogen ionisation 
front and the half-Irh point. Top left The H + ,H° and H2 density structure. Top right The C ionisation structure. Bottom left The O 
ionisation structure. Bottom right The Si ionisation structure. 



pie is the new grain model that has recently been devel- 
oped for CLOUDY (van Hoof et al. 2004). The major up- 
dates include a more flexible structure, allowing for grain 
sizes and species to be treated individually, rather than 
using a mean grain to be representative of an ensemble. 
This is important if the contribution of small grains to 
the SED is to be accurately predicted. Furthermore, these 
authors presented a more sophisticated treatment of the 
grain charge problem, inspired by the work of Weingart- 
ner & Draine (2001); the hybrid grain model improves on 
the previously implemented average grain charge model 
(Baldwin et al. 1991), allowing a more accurate calcu- 
lation of photoelectric yields. Figure 2 (van Hoof et al. 
2004) shows the effects of dust on the emission line spec- 
tra of typical PNe and H II regions models, obtained by 
including a dust component to the Meudon/Lexington 
benchmarks. The models in the plots show some dra- 
matic departures from the dust-free case; it should be 
noted, however, that these models extend to the PDR re- 
gion, while the effects on the photoionised region alone 
are expected to be smaller in magnitude. 



A new version (2.0) of the 3D Monte Carlo photoion- 
isation code MOCASSIN has recently been released (Er- 
colano et al. 2005), which includes a fully self-consistent 
treatment of the dust radiative transfer within the ionised 
region. Dust grains are allowed to compete with the ions 
in the gas for the absorption of UV and resonance emis- 
sion line radiation, with scattering by the grains being 
fully implemented in the isotropic approximation or via 
the application of a phase function. The microphysi- 
cal processes listed above, which further couple the two 
components, are also accounted for in the thermal and 
ionisation balance equations. The average grain charge 
model was preferred by these authors to a fully quantised 
or hybrid treatment in order to limit overheads, in light of 
the large uncertainties still existing in the dust data (see 
Weingartner & Draine 2001, Sec. 2.3). The MOCASSIN 
code, which can also be run in a pure-dust mode, has 
the advantage of being fully three-dimensional, hence al- 
lowing the scattering problem to be treated. Work is cur- 
rently in progress to allow polarisation maps to be pro- 
duced, which may yield better constraints on the geom- 
etry of the many sources for which this type of observa- 



tions are available. 



Expansion to Photon Dominated Regions 

Photon Dominated Regions (PDRs - also known as 
photo-dissociation regions) are associated with all pho- 
toionised nebulae, including PNe, that are optically thick 
to 13.6 eV radiation. When all photons with hv > 
13.6 eV have been absorbed H becomes predominantly 
neutral, in general marking the beginning of the PDR 
(Tielens & Hollenbach, 1985). Classically, the two re- 
gions have been treated as separate problems, both in 
terms of empirical analyses based on spectroscopic ob- 
servations and in terms of the development of numerical 
codes that would be able to model either of the two envi- 
ronments. In reality, a photoionised region and its PDR 
are closely coupled and, therefore, a self-consistent cal- 
culation is preferable. The radiation field impinging on 
the edge of a PDR is, in fact, a result of photon diffu- 
sion through the ionised region, and is naturally calcu- 
lated by a standard photoionisation code, whilst classical 
PDR codes generally treat Go 1 as a free parameter to be 
determined from the model. Moreover, some of the most 
important IR diagnostics used to characterise the gas in 
a PDR are also partially emitted in the ionised region. 
This is the case, for example, for all transitions of ions 
such as O , C + , Si + , having Ionisation Potentials (IPs) 
lower than 13.6 eV. Figure |2] shows a plot of the ionisa- 
tion structure of an H II region and its associated PDR, 
as recently calculated by Abel et al. (2005). 

Many PNe are known to be surrounded by extended 
PDRs, as is the case, for example, of NGC7027 (e.g. 
Graham et al. 1992, Latter et al. 2000) and NGC6302 
(e.g. Liu et al. 2001a). More recently, PDRs have also 
been detected around cometary knots, such as those of 
the Helix nebula (O'Dell, Henney & Ferland 2005) and 
the Ring nebula (Speck et al. 2003). 

The CLOUDY team are again leading the way, with 
a new version of their code (C05.07.06) recently re- 
leased, which allows self-consistent calculations of 
the spectrum, chemistry, and structure of photoionised 
clouds and their associated PDRs to be performed (Abel 
et al. 2005). Besides the improvements to the grain 
model mentioned in the previous section, a molecular 
network was also included (Abel et al. 2005) with 
approximately 1000 reactions involving 68 molecules. 
Recent advances in the treatment of H2 are described 
by Shaw et al. (2005). The chemical equilibrium 
was benchmarked at the Leiden 2004 PDR workshop 



Go is defined as the Far Ultra Violet (FUV) field between 6 and 
13.6 eV relative to the background interstellar radiation field (Habing, 
1968) 



( http://hera.phl. uni-koeln.de/^roellig/ Roellig et al. in 
prep.) and found to be in agreement with the results of 
other PDR codes. 

Having completed the inclusion of dust grains and 
subsequent benchmarking phase, the MOCASSIN code is 
also currently being expanded to incorporate the molec- 
ular reaction network in use by the UCL_PDR code (Bell 
et al. 2002, and references therein), which is amongst 
the codes recently benchmarked at the Leiden 2004 PDR 
workshop. The new version of the MOCASSIN code, 
which will be described in a forthcoming paper by Er- 
colano, Bell, Barlow & Viti., will allow one to make self- 
consistent models of ionised regions and PDRs of arbi- 
trary geometry and density distributions, including, for 
example, cometary knots and their shadow regions. 



Development of 3-D codes 

Available computer power has increased enormously 
since the development of the first generation of pho- 
toionisation codes in the late sixties. This has allowed 
the construction of more complex models, that, thanks 
also to the great advances in atomic physics, may include 
more ions, more frequency points, more lines and more 
atomic levels. Nevertheless, with a few exceptions, the 
fundamental assumption of spherical symmetry has al- 
ways been retained. A glance at spatially resolved images 
of H II regions and PNe, that are readily obtainable with 
modern instrumentation, immediately demonstrates that 
these object are rarely circular in projection. The effects 
of geometry and/or density distribution on the emerging 
emission line spectrum are in some cases dramatic, and 
an unrealistic model may result in the incorrect determi- 
nation of nebular properties (e.g. elemental abundances) 
and/or central star parameters, which may yield mislead- 
ing conclusions. 

The careful combination of a number of volume- 
weighted spherically symmetric models (composite 
models) has been used to treat some aspherical ge- 
ometries, including bipolar structures (e.g. Clegg et al. 
1987) or density/chemical enhancements (Pequignot 
2003). Morisset et al. (2005) have recently presented a 
pseudo-3D photoionisation code, based on this principle. 
The code consists of an IDL(RSI) suite of programs that 
make several calls to a ID code for a number of selected 
lines of sight through a 3D structure. The ID code 
currently used is the NEBU code (Morisset & Pequignot 
1996, Pequignot et al. 2001), but the technique can be 
readily applied to any other ID photoionisation code. 
Whilst this method has the advantage of being consid- 
erably faster than running a full 3D simulation, one of 
its major drawbacks is that the transfer of the diffuse 
component of the radiation field cannot, in general, be 



treated self-consistently. Ercolano et al. (2003b) found 
appreciable discrepancies between the results obtained 
by a self-consistent calculation of even a simple bi- 
conical geometry and those obtained by a composite 
model. 

To my knowledge the first two 3D photoionisation 
codes were those developed by Baessgen et al. (1990) 
and by Gruenwald et al. (1997). The latter is still in 
use and 3D models of aspherical PNe obtained with this 
code have recently been published (Monteiro et al. 2005, 
Monteiro et al. 2004). Both codes used analytical ap- 
proaches to the transfer of radiation; the first treats the 
diffuse field in the on-the-spot approximation, whilst the 
second employs an iterative technique for the determina- 
tion of the diffuse field. 

The fast development of parallel computing seen in 
last few year, has made the application of stochastic tech- 
niques to the radiation transport problem an attractive 
option. The first ID Monte Carlo photoionisation calcu- 
lations were those of Och et al. (1998). The first fully 
3D Monte Carlo photoionisation code, MOCASSIN (Er- 
colano et al. 2003a), uses a similar approach, but a more 
general version of the Monte Carlo estimators for the ra- 
diation field, developed by Lucy (1999). The description 
of the radiation field in terms of discrete quantities allows 
one to simulate the individual scattering, absorption and 
re-emission events that characterise a photon's trajectory 
as it diffuses through a cloud. This public code is fully 
parallelised and looks to exploit the increasing accessi- 
bility of large Beowulf clusters to allow the construction 
of realistic nebular models. As mentioned above, the in- 
clusion of grain physics has been one of the major recent 
updates of the code (Ercolano et al., 2005), in preparation 
for the expansion to PDR, planned for the near future. 

Wood, Mathis & Ercolano (2004) also presented a 
3D Monte Carlo photoionisation code, based on similar 
techniques as those described by Och et al. (1998) and 
Lucy (1999). This code is tailored for the modelling of 
Galactic H II regions and the percolation of ionising 
photons in diffuse ionised gas (Wood & Mathis, 2004) 



3-D photoionisation models of PNe 

In this section some recent results of 3D pure (no dust) 
photoionisation calculations recently obtained with the 
MOCASSIN code are summarised. 

Figure [3] (Neal, Ercolano & Sahai, in preparation) 
show the projected maps obtained from a model of the 
young PN MyCnl8, compared with the Hubble images 
presented by Sahai et al. (1999). The model images are 
shown to be in very good agreement with the obser- 
vations, indicating that the model 3D density distribu- 



TABLE 1. Model predicted emission lines 
for NGC 7009. Subset of model spectrum ob- 
tained by Goncalves, Ercolano et al., 2005a. 



Line Identification (A) No Slit Slit 



He I 5876 


14.1 


13.4 


He II 4686 


16.8 


22.0 


[O II] 3726,29 


5.38 


5.88 


[O III] 5007 


1197 


1172 


[N II] 6584 


5.87 


6.64 


[S II] 6731 


0.84 


0.99 


[S m] 6312 


1.47 


1.55 



tion for this object, previously presented by Dayal et al. 
(2000) and employed in our simulation is a realistic de- 
scription for this object. A well constrained density dis- 
tribution adds confidence to the model predictions of the 
other nebular and stellar physical properties. 

The effects of observations with a finite aperture on 
non-spherically symmetric objects should also be taken 
into account, their magnitude being larger the greater the 
departure from spherical symmetry of the observed ge- 
ometry. As an example, long-slit spectroscopic obser- 
vations of extended objects, such as H II regions and 
PNe, are routinely performed. Whilst line ratio measure- 
ments obtained from the summed spectra of a long-slit 
scanned across the whole 2D projection of a nebula may 
be considered to be representative of the integrated neb- 
ular spectrum, they are not always available. Geometry 
effects on long-slit spectra at a single PA were recently 
studied in the case of the PN NGC 7009 (Goncalves, 
Ercolano et al., 2005a, Goncalves et al., 2005b, these 
proceedings) and were found to be partially responsi- 
ble for the long-standing apparent overabundance of N 
in the FLIERs of this object (Balick et al. 1994, Hajian 
et al. 1997, Balick et al. 1998, and references therein). 
Table^lists a subsection of the slit-less (column 2) and 
slit (column 3) MOCASSIN model emission line spectra 
of the FLIERs of NGC 7009, showing an increase of ap- 
proximately 13% in the predicted flux of [N II] A6584. 
These results, discussed in more detail by Goncalves, 
Ercolano et al., 2005a, show that care should be taken 
when applying ICF methods or unrealistic photoionisa- 
tion models to non-spherically symmetric objects. In this 
light, even the classically accepted type i/ll classification 
scheme of PNe may become questionable. 

Chemical inhomogeneous models of the hydrogen- 
deficient polar knots of the PN Abel 30 were presented 
by Ercolano et al., 2003c; these models showed that it is 
in fact possible to explain the extreme ORL-CEL abun- 
dance discrepancy factors (ADFs) reported for this object 
by a simple bi-abundance model, consisting of a dense 
metal-rich (cold) core of ionised gas mixed with dust 
grains surrounded by a less dense envelope of hydrogen 



FIGURE 3. Left Panel: From Sahai et al. 1999. Narrow-band HST WFPC2 images of MyCn 18: (a) F658N ([N Il]16586; (b) 
F656N (Ha); (c) F502N ([O IIl]^.5007); (d) a continuum filter F547N. Right panel: From Neal et al. (in prep.). 2D projections of 
3D emissivity grids obtained from our best-fitting MOCASSIN model of MyCnl8: (a) [N Il]A6586; (b) H/3); (c) [O IIl]A5007; (d) 
[O i]A6300. 



deficient gas with less extreme metal abundances. In this 
model, ORLs of heavy elements result are shown to be 
emitted almost exclusively by the cold core, CELs from 
the hotter surrounding envelope, whilst He recombina- 
tion lines are emitted by both components. 

In the case of some H II regions for which only mild 
ADFs are observed, the discrepancy may be resolved 
by considering temperature fluctuations in the t 2 pre- 
scription of Peimbert (1967). Whilst this interpretation 
has been successful in explaining the observations for 
some objects with mild ADFs, the nature and cause of 
the supposed temperature fluctuations remain uncertain. 
Even models that include density fluctuations yield only 
small deviations from an otherwise near homogeneous 
temperature structure of the emitting regions. One pos- 
sibility worth investigating is whether larger fluctuations 
may be produced with a model including multiple ion- 
ising stars. A number of 3D MOCASSIN models are be- 
ing constructed which use several geometries and ionis- 
ing source distributions (Ercolano, Stasihska & Barlow, 
in preparation), which may help to clarify the situation 
from a theoretical point of view. 

CONCLUSION: THE FUTURE 

It has been the aim of this review to identify and dis- 
cuss in broad terms the most recent advances in pho- 
toionisation codes. Several large scale codes exist and are 
being maintained by individual working groups, coming 



together regularly to ensure the reliability of the results 
via the performance of a set of benchmark problems. 
Modern codes aim at the construction of realistic sim- 
ulations in order to efficiently extract information from 
a wealth of new spectroscopic observations available at 
large wavelength ranges and for increasingly more dis- 
tant and faint objects. 

The coupling of photoionisation/PDR and hydrody- 
namic RT codes is only a step away. This is perhaps 
one of the most exciting prospects, as the new simula- 
tions will undoubtedly yield new insights into the physics 
of many astronomical phenomena. Already some pho- 
toionisation codes allow for time-dependent effects to 
be taken into account. In the near future it will be pos- 
sible to perform a full photoionisation/PDR calculation 
at each time-step in a large hydrodynamic simulation, to 
follow, for example, the evolution of the time-dependent 
non-equilibrium spectrum of a dynamically evolving gas 
cloud, from the time it starts being illuminated. Both 
physics and numerical techniques are mature, it is only 
a matter of waiting for the computer power required to 
become available. 
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